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Abstract 

We consider the computation of D and B fragmentation functions in e + e~ 
annihilation. We compare the results of fitting present data using the next- 
to-leading- logarithmic resummed approach, versus the 0{a 2 s ) fixed-order cal- 
culation, including also mass-suppressed effects. We also propose a method for 
merging the fixed-order calculation with the resummed approach. 
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1 Introduction 

A theoretical framework for the study of heavy-quark fragmentation functions 
(HQFF) has been available for a long time [|l[], and several phenomenological analy- 
sis based upon this formalism have appeared in the literature |2|, |3|, |J. In Ref. 0, 
next-to-leading fits to the charm-momentum spectrum at ARGUS were performed 
and used to predict the bottom spectrum from Z decays. From this analysis, the 
non-perturbative part of the fragmentation function predicted for the bottom quark 
turned out to be quite hard: in fact, much harder than predicted in Monte Carlo 
models using the standard Peterson || parametrization. More precise data || [?]. § 
have become available since then, giving an indication of a hard bottom fragmenta- 
tion function. In Ref. || , a study of the charm fragmentation function was performed, 
using a parametrization of the non-perturbative effects based upon the Peterson frag- 
mentation function, instead of the form adopted in Ref. ||. A quantitative result on 
the value of the e parameter was obtained there, definitely showing that e is much 
smaller in next-to-leading-log (NLL) fits rather than in the leading-log (LL) ones. 

The interest in further refining our understanding of the fragmentation functions 
for heavy quark stems mainly from the possibility of using them to improve our un- 
derstanding of heavy-flavour hadroproduction and photoproduction. In Ref. ||, a 
formalism for the computation of heavy-flavour production, that merges the high- 
transverse-momentum approach (i.e. the fragmentation-function approach) and the 
fixed-order one, was developed. It was found there that mass-suppressed effects are 
quite large even at moderately large transverse momenta. This work is a first step 
towards a sound application of the fragmentation-function formalism in the moderate- 
transverse-momentum range. It should, however, be complemented by similar calcu- 
lations in the context of e + e~ annihilation, since this is the place where the impact 
of non-perturbative effects is studied. This is in fact the purpose of the present work: 
to use all the available knowledge on heavy-flavour production in e + e~ annihilation 
in order to reach an assessment of the size of non-perturbative effects. 

There are theoretical approaches to the fragmentation-function calculation that 
rely upon heavy-quark effective theory in order to study more systematically non- 
perturbative effects fpL0| . In the present work, however, we want to establish a con- 
nection with the most commonly used parametrization, and thus we will use the 
Peterson form throughout. 
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Fixed-order (FO) calculations of the O (a 2 ) differential cross section for heavy- 
quark production in e + e~ annihilation JD], [12|, |13|| do exist today, and some appli- 
cations in the context of HQFF [|14| , pL5 ] have already appeared. In practice, the 



fixed-order calculation should be more reliable than the HQFF approach for small 
annihilation energies. It is interesting, therefore, to consider an approach in which 
the fixed-order and the HQFF calculations are merged, in the spirit of the work of 
Ref. 0, without overcounting. In the following, we will thus review the HQFF and 
the fixed-order-calculation results on fragmentation functions. We will then define a 
merged approach, in which the fixed-order result is supplemented with leading and 
next-to- leading logarithmically-enhanced contributions at all orders in a s . We will 
consider both charm- and bottom-production data, which we fit using a Peterson 
parametrization of the non-perturbative contribution. 

The paper is organized as follows. In Section [| we fix our notation, and describe 
the aim of our improved approach. In Section [| we describe the procedure we use 
in order to incorporate a parametrization of the non-perturbative effects in our for- 
malism. This procedure is to some extent arbitrary. It is however important that 
the same procedure is used throughout our calculation. In Section [| we describe the 
fixed-order calculation of the fragmentation function. Some subtleties, related to its 
normalization, are discussed in detail. Section [5] is dedicated to the computation of 
the resummed cross section truncated to order a 2 , at the NLL level (TNLL). This 
calculation is almost equivalent to neglecting all mass-suppressed terms (i.e. terms 
that vanish like powers of m 2 /E 2 ) in the FO calculation, except that terms of order 
a 2 , without any power of log E 2 /m 2 , are not included here. In this section we spell 
out the definition of the improved approach, and clarify its practical implementation. 
In Section [] we compare the various approaches, and in Section [7| we describe our fits 
to various data sets. Finally, in Section ||, we give our conclusions. 

2 Theoretical framework 

We consider the inclusive production of a heavy quark Q of mass m 

e + e-^Z/ 1 (q)^Q(p)+X , (2.1) 



-4- 



where q and p are the four-momenta of the intermediate boson and of the final quark. 
Defining; the scaled energy of the final heavy quark 

2p- q 



X E 



Q 



(2.2) 



and introducing the centre-of-mass energy E = y/cp, we have the physical constraint 

y/p<x E <l, (2.3) 

where 

P=*£- (2-4) 



E 2 

At times, we will instead use the normalized momentum fraction, defined as 



x r 



p 



0<x p <l. 



VT-P' ----- (2 ' 5) 
The inclusive cross section for the production of a heavy quark can be written as a 
perturbative expansion in a s 
da 



x p , E,m) = ^2 a i x p) E, m, ji) a™ (//) 



"P n=0 

where \x is the renormalization scale, and 



a s (fi) 



2tt 



(2.6) 



(2.7) 



If fx w E w to, the truncation of Eq. ( |2.6|) at some fixed order in the coupling con- 
stant can be used to approximate the cross section. An O (oil) fixed-order calculation 
for the process ( |2.1| ) is available [ HI O, O, so that we can compute the coefficients 
of Eq. (|2~6|) at the (9 (a^) level. We thus define the fixed-order result as 

da 



(x p ,E,m) 

UjJU <n 



a {0) (x p ,E,m) + a w (x p ,E,m)a s (E) + a (2) (x p ,E,m)a 2 s (E) , (2.8) 



(2), 



FO 



where we have taken \i — E for ease of notation. 

If E ^> m, large logarithms of the form log (E 2 /m 2 ) appear in the differential 
cross section ( |2.6j ) to all orders in the perturbative expansion. In this limit, if we 
disregard all power-suppressed terms of the form m 2 /E 2 , the inclusive cross section 
can be organized in the expansion 



n=0 



(»), 



E 



x) <x 



(E)\og—\ +a s (E)^l 



x 



m 



a s (E) log ^ 



n=0 



+ a 2 s (E)Y / S {n) (x) la s (E)\og 



n=0 



E 2 



+ ... + 



m 
E 2 



(2.9) 
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where x stands now for either or x p , since the two variables differ by power- 
suppressed effects. 

We define the leading-logarithmic (LL) approximation as 



da ^ 

dx 



x, E, m) 



£/? ( 



n). 



X) 



a s {E)\og 



E 2 



LL n=0 \ 

and the next-to-leading-logarithmic (NLL) one as 



77T 



(2.10) 



da ^ 

dx 



x, E, m 



E 



NLL n=0 



J2P in) (x) a s (£)log — +a s (E)Y / l {n) (x) [a s (E)\og 



E 2 



77T 



12.11) 



The 5^ coefficients define the NNLL terms, that of now, not known. 

The expansion of Eq. (|2.11|) up to order a 2 is given by 



da 

dx 



'x, E, m) 



P i0 \x) + a s {E) I 7 w (x) + (3 w {x) lo 



(o), 



NLL 



E 2 



+ a 2 s (E)[^\x) 



log^ \ + fl 2 \x) log 2 &)+ (a£) , (2.12) 



and does not coincide with the massless limit of the FO calculation. A term of order 
a 2 , not accompanied by logarithmic factors, may in fact survive in the massless limit 
of the FO result. In the HQFF approach, this is a NNLL effect, and therefore it is 
not included at the NLL level. We will refer to this term, in the following, as the 
NNLL a 2 term. 

It is now clear how to obtain an improved formula, which contains all the informa- 
tion present in the FO approach, as well as in the HQFF approach. Using Eqs. (|2. 
and (|2.11|) , we write the improved cross section as 



da 

—{x p ,E,m) 

UiJb ft 



J2a {i) (x p ,E,m) al(E) + J2 P (n) 



x 



' E 2S 
a s (E) log — 



imp 



i=0 



n=3 



(n). 



X) 



f E 2 N 

a s (E) log — 



n=2 



(2.13) 



where the LL and NLL sums now start from n = 3 and n = 2 respectively, in order 
to avoid double counting. Formula fl2.13|) includes exactly all terms up to the order 
a 2 (including mass effects), and all terms of the form (a s {E)\og{E 2 /m 2 )) n and 
a s (E) (a s (E) log (E 2 /m 2 )^J , so that it is also correct at NLL level for E ^> m. It 
can be viewed as an interpolating formula. For moderate energies, it is accurate to 
the order a 2 , while for very large energies it is accurate at the NLL level. 
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3 Non-perturbative effects 

In the HQFF formalism, the inclusive heavy-flavour cross section is given by the 
formula 

^-(x,E,m) = ^ (x, E, fj, F ) <g> L>i (x, fj, F , m) 

i 

~ l dz da. 



^ f l dz dai -fx \ 

2_j] — {z, E, n F ) Di I — , [j, F , m J , (3.1; 



where d&i(x, E, fi F )/dx are the MS-subtracted partonic cross sections for producing 
the parton i, and Di(x, ^i F ,m) are the MS fragmentation functions for parton i to 
evolve into the heavy quark Q. The factorization scale \i F must be taken of order 
E in order to avoid the appearance of large logarithms of E/fip in the partonic 
cross section. The explicit expressions for the partonic cross sections and for the 
fragmentation functions at NLO can be found in Refs. [1], [16|]. 



The weak point of formula ( |3.1|) comes from the initial condition for the evolution 
of the fragmentation function, which is computed as a power expansion in terms of 
a s {m). In particular, irreducible, non-perturbative uncertainties of order A QCD /m 
are present. We assume that all these effects are described by a non-perturbative 
fragmentation function D$ P , that takes into account all low-energy effects, including 
the process of the heavy quark turning into a heavy-flavoured hadron. The full 
resummed cross section, including non-perturbative corrections, is then written as 

-^-(x,E,m) = J2^(x 7 E,fi F ) (2) Di(x,n F ,m) <g> D$ P (x) . (3.2) 
ax i dx 

We will parametrize the non-perturbative part of the fragmentation function with the 
Peterson form 

Dl v [x) = P(x, e) = N r * (1 ~* )2 l2 , (3.3) 

[(1 — x) z + xe\ 

where the normalization factor N determines the fraction of the hadron of type H in 
the final state. Summing over all hadron types, we have the condition 

£ [ 1 dxD^ P (x) = l. (3.4) 
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4 The fixed-order approach 

We define our fixed-order cross section, supplemented with non-perturbative frag- 
mentation effects, by the convolution of the perturbative cross section fl2.8|) with the 
non-perturbative fragmentation function ( |3.3|) 
da H 



(x p , E, m) 

( I ■( 11 



f 1 , i da 

J dydZp —{z p ,E,m j 



2 



P(y,e)6(x p - yZp) 



FO 



= / dydz p J2 a {l Xzp,E,m)al{E) P(y,e) 8(x p - yz p ) .(4.1) 

Jo i=0 

It is assumed, in this formula, that non-perturbative effects degrade the momentum, 
rather than the energy of the heavy quark. This is to avoid unphysical results when 
the momentum is smaller than the mass. It should be made clear, however, that for 
small momenta this formula should be seen merely as a model for non-perturbative 
effects. 

The coefficients of the perturbative expansion a^(z p , E,m) are in general distri- 
butions in z p with singularities at z p — 1. In order to deal with these singularities, 
we transform Eq. (|4.1|) as follows 

da H 



r 1 ■ 

= / dydz P Y,a M {zp,E,m)a l s {E) P(y,e) {5{x p -yzp) - 5{x p -y)} 
Jo i=0 

+ / dydz p J2a {l) {zp,E,m)al{E) P{y,e)5{x p -y) 
Jo i=0 

= P(xp,e)j2d^(E,m)al(E) 

i=0 

r 1 2 

+ / dydz p ^2a {l) {zp,E,m)al{E) P{y,e) {5{xp-yzp)-5{x p -y)},{4.2) 



where 



i=0 

a^{E,m)= dzpa (i) (zp,E,m) . (4.3) 
Jo 

The term in the last integral of Eq. ( |4.2[ ), in fact, does not contribute, since it is 
proportional to 5(1 — z p ). We thus find 
da H 



P(xp, e) o"; nc 



FO 



r 1 2 

+ / dydzpJ2a {i) {zp,E,m)al{E) P{y,e) {5{xp-yz p ) -5{x p -y)} ,{AA) 
Jo i=i 
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where 

a inc = a {0) (E,m) + a ( - 1) (E,m)a s (E) + a ( - 2) (E,m)a 2 s (E) (4.5) 

is the inclusive heavy-quark cross section. 

It would be natural to normalize the cross section in terms of the total cross section 
for heavy-flavoured events. However, for practical purposes, the normalization is 
immaterial, because of the uncertainty in the specific heavy-flavoured hadron fraction, 
and it must be fitted to the data. It is important, however, that the normalization 
factor has a finite massless limit, since the same normalization has to be applied to 
the resummed cross section. In the present case, neither the total heavy-flavour cross 
section nor the inclusive cross section are finite in the massless limit. This problem 
arises because of mass singularities in the coefficient a^(E,m). We then define 

a^(E,m) = a$(E,m) + 2a$(E,m)+a%\E,m) , (4.6) 

where: 

- the lQ-term consists of all graphs in which the primary interaction vertex is 
always attached to the heavy flavour, and furthermore there is only one heavy- 
flavour pair in the final state. 

- The 2Q-term consists of all graphs in which the primary interaction vertex is 
always attached to the heavy flavour, and there is a secondary heavy-quark pair 
in the final state coming from the gluon-splitting mechanism. There is a factor 
of 2 in front of this contribution, since, in the inclusive cross section, both the 
primary heavy quark and the secondary one can be detected. 

- The .R-term (where R stands for "rest" ) contains all other contributions. These 
include terms in which a heavy-flavour pair is produced via gluon splitting in 
a process initiated by light quarks, plus interference terms in which a heavy 
quark (antiquark), produced via gluon splitting in an amplitude, interferes with 
a heavy quark (antiquark), produced directly at the Z/j vertex. 

We define a "normalization" total cross section as 



<t„ = a w (E,m) + a w {E,m)a s {E) + [a$(E, m) + a^E, m)\ a' 2 s (E) . (4.7) 

This cross section is finite in the massless limit. In fact, the mass singularities, present 
in the a 2 Q(E, m) term, cancel against the virtual graphs in o\q(E, m) having a gluon 
self-energy correction consisting of a heavy-flavour loop. 
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The analytic expressions of (z p , E, m) and a' 1 ' (z p , E, m), and of (E, m) and 
a^(E, m), are known (see Refs. JIB], |T7|). The a^(z p , -E, m) term has been computed 
by the authors and is available as a numerical FORTRAN program. 

The advantage of writing Eq. fl4.1|) in the form of Eq. (fOJ) is that the singularities 
in the two-jet limit are regularized, and the two-body virtual terms (i.e. Z/7 — > Q+Q), 
that are proportional to 5(1 — z p ), give zero contribution. Thus, the two-loop virtual 

corrections, which were not computed in Ref. ||11|| , do not contribute there. On the 

(2) 

other hand, they contribute to a\Q(E,m). They do disappear from our formula, 
however, if we normalize it to the cr n cross section, consistently dropping higher order 
terms. We obtain 



1 



da 



H 



(x p , E, m ) 



FO 



P(x p , e) n R + _ 1 f dy dz p \ a {1 \z p , E, m) a s 

a^ u >{E, m) Jo \ 



+ 



a 



(i) 



(E,m) 
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a {1 \z p ,E,m) 
a (,J >{E, m) 



a 2 . 



x P(y, e) {S(x p - yz p ) -5(x p -y)} 



where 



1 + 



a™{E,m) + a%\E,m) ^ 
a(°)(E,m) 



o? s {E) . 



(4.9) 



(2) 

Although the d R (E, m) term contains all sort of interference contributions, it is dom- 
inated by terms in which a primary light-quark pair generates a gluon that decays into 
a heavy-quark pair. This contribution is singular in the massless limit. An analogous 
singularity is present in the d^Q(E,m) term. Thus, the dominant corrections to tir 
arise from secondary heavy-quark production via gluon splitting. 

Performing the integration in y we can write 

1 da H 



[x p , E, m 




dz p b (i) (z p ,E, m) 



P(xp/z p ,e) 



- P(x p ,e) 



P dzpb^(z p , E,m) 



(4.10) 



where we have defined 
b^(z p ,E,m) = 



a^'(z p , E, m) 
aW(E,m) 
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b®(z p ,E,m) 



(z p , E, m) - a (1) (zp, E, m) 



a(°)(E,m) 



■ (4.H) 



For future use, we introduce here the hadronic differential cross section expressed 
in term of the energy fraction. From Eq. ( |2.5|) , we can write 

1 da H 1 da H dx v 1 da H Xe 1 „„ s 

= L = (4-12) 

a a dx e c" n dx p dx e o" n dx p x p 1 — p 

5 NLL higher-order effects 

According to Eq. (|2.13| ), we must add to the FO result all terms of order oP s 
and higher of the full NLL resummed result. The resummed result is obtained nu- 
merically by solving the Altarelli-Parisi evolution equations for the fragmentation 
functions, with the given 0(a s ) initial conditions, convoluted with the appropriate 
short-distance cross sections. We must, however, subtract, from this numerical result, 
the terms up to second order of its expansion in powers of a s . These terms can be 
obtained with the same procedure used in Ref . ]14| . We recall here the principal steps 
that give us the O (a 2 s ) terms of the resummed cross section. 

We introduce the following notation for the Mellin transform of a generic function 
/(*): 

/(AO = / dxx N ~ l f{x) . (5.1) 
Jo 

We adopt the convention that, when iV appears, instead of x, as the argument of 
a function, we are actually referring to the Mellin transform of the function. This 
notation is somewhat improper, but it should not generate confusion in the following, 
since we will work only with Mellin transforms. 

The Mellin transform of the factorization theorem ( p.l| ) is given by 

a{N, E,m) = J2 &i(N, E, /i) A (A, y., m) , (5.2) 

% 

where 

a(N, E, to) = f 1 dx x N ~ l -^(x, E, to) , (5.3) 
jo dx 

and a similar one for <7j(A, E, jj). The Mellin transform of the Altarelli-Parisi evolu- 
tion equations is 

dDt l!ll% m) = £«.(A0 [ P vW + a.M + O (cS)] Dj{N,y,m) . (5.4) 
a log /i ■ 
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We need an expression for a(N, E, m) valid at the second order in a s . Thus, we solve 
Eq. (5^4), with initial condition at /i = /io, with an O (a 2 ) accuracy. This is easily 



done by rewriting Eq. (|5.4j) as an integral equation 
Di(N,(i,m) = Di(N,iio,m) 

+ E rdlogfi' 2 a a (ji) [4 0) (iV) + P£\N) a s (//)j D^N.^m) . (5.5) 



The terms proportional to ot 2 can be evaluated at any scale (// or /i ), the difference 
being of order oP s . Factors involving a single power of a s can instead be expressed in 
terms of a s (fio) using the renormalization group equation 

a s (/i') = a s (fi ) - 2 -k b a 2 (/i ) log — + O (a 3 s (^ )) 

Mo 

UC A -4n f T F 

b o = T7T- ~ > (5-6) 

LA 71 

with nf the number of flavours, including the heavy one. Equation ( |5.5| ) then becomes 
A.(iV,/i,m) = A(iV,// ,m)+^ I* d\ogti' 2 a a {to)P§\N)D j {N, l i\m) 

+ E "'(Mo) P}P(N) Dj(N , no, m) log ^ 

3 ^0 

-2nb E^(Mo) Pij\N) 4(iV,// , m)\ log 2 ^ . (5.7) 
i ^ Mo 

Now we need to express Dj (N, //, m) on the right-hand side of the above equation 
as a function of the initial condition, with an accuracy of order a s . This is simply 
done by iterating the above equation once, keeping only the first two terms on the 
right-hand side. Our final result is then 

2 

Di(N, fi, m) = Di{N, /i , m) + E a s (fi )P^(N) Dj(N, // , m) log ^ 

j Mo 

+ E «^o)4 O) (A0 Pg } (iV) 4-(iV, aio, m) \ log 2 ^ 
fcj z Mo 

+ J2^o)P}p(N)D j (N, f , ,m) log ^ 
j Mo 

- vr 6 E a^^V) 4-(^> Mo, m) log 2 ^ • (5-8) 

Mo 
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Using the initial conditions (see Ref. 0) 

D Q (N,ij ,m) = l + a s (/io)rfJ 3 1) (A^,/i ,m) + C(a^(/io)) 

b g (N,fi ,m) = a 8 (jio)SJ-\N,fj^,m) + 0(o^(jio)) , 

(all the other components being of order a^), and using Eq. ( |5.6|) to express a s (/x ) 
in terms of a s (//), Eq. (|5.8|) becomes, with the required accuracy, 



(5.9) 



Di{N,n,m) = 5 iQ + a s (p) df\N^ Q) m) + 2vr6 a^(/i) d\ L> (N, /j, ,m) log^ 

+ a s (f,)P$(N) log ^ + E &Mlf{N) A*o, m) log 4 

+ E * 2 M4 o) (A0 p!3(n) \ lQ g 2 4 + &MP$W log 4 



/'0 



+ 7r6 a s 2 (/i)P J g ) (iV) lo. 



5 2 
00 



Using the notation of Ref. 16 1, we have, for the partonic cross sections, 



°Q = °q = a o,Q 1 + a s C F c Q + O (afj 
bq = a q = <To, g [1 + (a s )} 



cr g = cr 



a s C F Cg + O (a 



where Q is the heavy quark and q is a light one, and 

a 0,Q 



_(«) , _(«) 



(v) . (a) 
a 0,g + <?0,q 



/'6 



(5.10) 



(5.11) 



°"0 = E °"0J = C"0,Q + E% 
CQ = C T ,Q + C L ,Q 

Cg = C T ,g + C L ,g, (5.12) 

where v, a, T and L refer to the vector, axial, transverse and longitudinal contribu- 
tions, respectively. 

In addition, we define the Mellin moments 



a§ g (N,E, l i) = C F c Q/g (N,E, t i) 



(5.13) 



13 



and introduce the total "normalization" O (a s ) cross section for the production of a 
heavy quark 

r 3 

°n = oo,Q l + -C F a s + 0(af) . (5.14) 

Observe that now, at NLL accuracy, we can neglect the O (a 2 ) terms, since they do 
not contain any logarithmic enhancement. This would not be the case for the total 
heavy-quark cross section, that in fact is divergent at order O (a 2 ), in the massless 
limit approximation. 



From Eq. ( |5.2|) , we can write the NLL cross section as 

a(N) = a Q D Q + aQDQ + Y / ^ q D q + Yl °H A + A? <^ > 



(5.15) 



where we have dropped the energy and the mass dependence, for ease of notation. 
We can now obtain the truncated O (a 2 ) NLL (TNLL from now on) normalized 



cross section as 
1 



-a{N) 



TNLL 



1 + a.(ji) d%\N) + aMPj&iN) log^ + <*M U$\N) - \c 



+ ^)P§kN) ( W{N) - ~-C F ) log ^ + 2 tt b of (/i) d%(N) log ^ 



;(!)/ 



t'-o 

2 



/(I)/ 



/'0 



/'o 



+ T,^)PQhN)d^(N) log^ 



+ E ^(A*)PS W P fe ( g(iV) 1 log 2 IL + afo)P« (JV) log £ 



+ nb a'MPQQ(N) log 



Mo 



+ &MP%W df\N) log g + a^)Pff g (N) P$>(N) \ log 2 ^ 



+ « 2 (/i)P«(iV) log _ j2 



Mo. 



+ ^E 5 -o, 9 « 2 (/i)<|pW(iV)4 1 )(iV) log^ 



a 0,Q 



/'o 



+ ^ (AO P$ (N) ~ log 2 £ + P# (iV) log ^ 



2 

/'0 
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+ pg)(N)dU(N) log 



EL 

f4 



Mo 



Mo 



+ -J-^af\N)a 2 s ^)P^(N) log ^ , 
^o.q ' Mo 



(5.16) 



where we have taken \i = E and /i = m, and we have introduced the notation 



h Q/M = ^ /a {N,E^)\, =E , d% /g (N) = d^ g (N^ ,m)\, 0=m 

The lowest order splitting functions are given by 



;(!) 



(5.17) 



iV(iV 



2S r i(iV) 



P%{N) = P£ j (N) = pg\N) = C F 



Qg 



2 + N + N 2 
N(N 2 - 1) 



where 



and 



2 + JV + N 2 
N(N + 1)(N + 2) 

Si(N)=<0b(#+l)-Vo(l) 
d 



log Hz) . 



In addition we have 



P$(N)=P$(N) = P q , q (N) 



(5.18) 

(5.19) 
(5.20) 
(5.21) 



and the splitting functions Pqq(N), P^(N) and P q > q (N), together with the ini- 
tial condition for the fragmentation functions dg (N) and d g ^(N), can be found in 
Refs. |, [n|. 

The differential cross section, as a function of x, can now be obtained by an inverse 
Mellin transform. We define the hadronic cross section, including non-perturbative 
effects 

1 do H 



<7 n dx 



(x, E, m 



TNLL 



1 fc+ioo 1 

2^ y ■ dN v a{N) 

LiK L J C — ZOO Or> 



P(N,e)x 



-N 



(5.22) 



TNLL 



where P(N, e) is the Mellin transform of the Peterson function. We can now isolate 
the higher-order effects (HOE) needed in the expression of the improved cross section 
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of Eq. (gl3|) . They can be computed as a difference between the NLL and the TNLL 
cross section 



00 / E 2 



(HOE) = E/5 (n) (^) «-(^)log 



n=3 



x / 



ra=2 



da 



H 



dx 



(x, E, m) 



NLL 



da H , „ s 



TNLL 



(5.23) 



where 



da 



H 



dx 



(x, E, m) 



NLL 



^(x,E,m) 
dx 



P(x,e) . 



NLL 



(5.24) 



We can then summarize our improved cross section as 



1 da H 



1 da H 



imp 



+ 



1 da H 



FO 



<7 n <ix 



1 d<j H 

NLL °" n ^ X 



TNLL 



(5.25) 



or a similar one with x p replaced by xe- 



6 Comparison of the various approaches 

In this section we perform a comparison of the different approaches presented 
so far. We confine ourselves to D and B meson production at E = 10.6 GeV and 
E = 91.2 GeV, which are relevant to the data sets that we will fit in the forthcoming 
section, where we analyze the experimental results obtained by ARGUS and OPAL 
for D mesons, and by ALEPH for B mesons. We fix the charm and bottom mass to 
1.5 and 5 GeV, respectively. For D meson production we take €d = 0.035, while for 
B mesons we use €b = {rn 2 c /m^) ~ 0.0035. The renormalization and factorization 
scales have been settled equal to the total energy E. Furthermore, we have fixed 
Aqqd = 200 MeV, so that af\M z ) = 0.116. 

In the fixed-order calculation, we have taken three light flavours, for the energy of 
10.6 GeV, and four light flavours for E = 91.2 GeV. According to the renormalization 



scheme []I8[ we used in our FO calculation pT[ , this implies that the strong coupling 
constant runs with (rif — 1) flavours, where rif is the total number of flavours, includ- 
ing the massive one. For this reason, in order to compare the FO results with the 
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resummed and the truncated ones (where the the heavy flavour is treated as a light 
one), we have to change the strong coupling constant, according to 



a 



(6.1) 



3 " ~ m c 

where we have used the renormalization group equation (|5.6|) and the matching con- 
dition 

'^(m). (6.2) 



(n/-l) 



This implies that the FO differential cross section, computed with a s 

-- a<°> + a« af^iE) + a® a 2 s (E) , (6.3) 



da 



FO 



acquires a contribution, proportional to the aW term, once written in terms of ai 1 ^ 



da 

dec t) 



a (o) +a (D «("/)(£) + 



FO 



iV- 2 -T F aV log^ 
3 m z 



ai(E) , (6.4) 



where we have dropped all the arguments of the coefficients, for ease of notation. 
In Figs. 0, Q and |3], we plot the following quantities: 

- the dot-dashed line is the differential cross section at order a s . It corresponds 
to the sum of the terms up to the order a s in Eq. ( |4.10|) ; 



the dotted line represents the cross section at order a s without power-suppressed 
mass effects. We see that such terms give a noticeable contribution only at 
ARGUS, while they are completely negligible at OPAL and ALEPH; 



the dashed line represents the LL resummed cross section of Eq. (|2.10|) . 



the solid line is the improved differential cross section, obtained by merging the 
LL resummed and the a s -order massive one. At OPAL and ALEPH the LL and 
the improved LL curves are almost identical. From this we infer that O (a s ) 
mass terms are small at ARGUS energies, and completely negligible both for b 
and for c at LEP energies. 



Notice that for very large and very small x, the distributions we computed may 
become negative. This is an indication of the failure of the perturbative expansion, 
due to the presence of large terms proportional to powers of log(x) and log(l — x). 
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These terms have not been resummed in our approach. We will discuss in more detail 
this problem in the next section. 

Figures [|, |5] and [| are similar to Figs. [I], |] and |3|, but they include next-to-leading 
effects. 

Thus: 



the dot-dashed line is the differential cross section of Eq. (|4.10| ) at order a 2 . All 
mass effects and the exact a 2 term are present; 



the dotted line represents the TNLL cross section of Eq. ( |5.22 ). It differs from 



the previous curve for the absence of mass power-suppressed effects and of NNLL 
a 2 contributions (the terms not accompanied by large logarithms). We can see 
that, for the OPAL and ALEPH curves, where m 2 /E 2 terms are much smaller 
than the ARGUS ones, the effect of the NNLL a 2 term is quite small, except for 
the small- a; region in the OPAL cross section, where the splitting mechanism of 
a gluon, coming from a primary light-quark pair, gives a sizable contribution; 



- the dashed line is the NLL resummed cross section of Eq. fl5.24Q , normalized to 
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X! 

n 

b 



2.0 



1.5 — 



1.0 — 



0.5 



0.0 




OPAL with e = 0.035 
Solid: LL improved 
Dashes: LL 

Dot-dashes: oc s fixed order 
Dots: a s fixed order with 

no mass-power terms \ 
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0.4 



0.6 



0.8 



Figure 2: Fragmentation function l/a to t da D /dxE for OPAL. The value of 
the Peterson e parameter is fixed at 0.035. The dashed curve is almost hidden 
by the solid one. 



cr n of Eq. (ETI3D 



- the solid line is the improved cross section of Eq. fl5.25| ), that takes account of 
the NLL resummed effects and of the mass and NNLL a 2 s terms. 

We notice that, for the next-to-leading curves, the resummed cross section tends to 
be softer than the fixed-order one, so that it will need a harder non-perturbative 
fragmentation function (corresponding to smaller values of e) in order to fit the data. 
Furthermore, this effect is much more pronounced at LEP energies, as one can expect. 

Mass effects seem to be very small in this context. In general, we see that they 
harden the fragmentation function. We thus expect that they will lead to larger 
values of e when fitting the data. 

7 Fit to the experimental data 

We present now some fits to the experimental data in order to extract the non- 
perturbative part of the fragmentation functions. 
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We consider the following data sets 



D* + meson at ARGUS 111, at E = 10.6 GeV 



- D* meson at OPAL |0fl, at E = 91.2 GeV, 

- B meson at ALEPH §, at E = 91.2 GeV. 

Besides the LL and NLL fits (similar to that ones made in Ref. [H]), we present new 
fits with our improved cross section. Furthermore, we present fits where the initial 
conditions ( |5.9| ) for the fragmentation functions are taken in an exponentiated form 
(called "NLL expon"). In iV-space, the exponentiated initial conditions read 



D Q (N, ii Q ,m) 
D g (N, ii ,m) 



exp 
exp 



a s (fi ) dft (N, no,m) 
a s (/j, ) df\N, Ho,m) 



1 . 



(7-1) 



These initial conditions are equivalent, from the point of view of NLL resummation, 
to the ones in Eq. ( |5.9|) . They are introduced here only to enhance higher order effects 
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in the large-x region. In fact, these terms partially account for Sudakov behaviors 
in the endpoint region |], [Hp. We will see that their effect is quite dramatic, and, 
in particular, that they remedy to the problem of negative cross sections at large x. 
They are shown here just for the purpose of illustrating how higher order perturbative 
terms may get rid of this problem. 

The negative values in the small-x region are related to the multiplicity prob- 
lem We will not try to remedy to them in this context, since they are not very 
important in the experimental configurations we consider. 

We have fitted the data by \ 2 minimization. With this procedure we have fitted 
both the value of e and the normalization, which was allowed to float. We have kept 
Aqqd fixed to 200 MeV. 

The results of the fits are displayed in Tables [l] and |2|. The corresponding curves, 
together with the data, are shown in Figs. 0-^. 

The full improved resummed result of Eq. ( |5.25| ) has been used here. For compar- 
ison, we have also plotted the NLL curves computed at the same value of e. 
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e (x 2 /dof) 


a s fixed order 


LL 


LL improved 


ARGUS D 


0.058 (0.852) 


0.053 (2.033) 


0.054 (2.194) 


OPAL D 


0.078 w (0.706) 


0.048 (1.008) 


0.048 (1.008) 


ALEPH S 


0.0069 (4.607) 


0.0061 (0.137) 


0.0064 (0.137) 



Table 1: Results of the fit of the non-perturbative e parameter for the Peterson 
fragmentation function. The value of x 2 /dof is given in parenthesis. The range of the 
fit is indicated in Figs. P|-IT with small crosses. '*'We have excluded the first three 
experimental points for the a s fixed-order fit to the OPAL data. 



e (xVdof) 


a 2 fixed order 


NLL 


NLL improved 


NLL expon 


ARGUS D 
OPAL D 
ALEPH B 


0.035 (0.855) 
0.040«(0.769) 
0.0033 (2.756) 


0.018 (1.234) 
0.016 (1.122) 
0.0016 (0.441) 


0.022 (1.210) 
0.019 (1.066) 
0.0023 (0.635) 


0.0032 (1.493) 
0.0042 (1.152) 
0.0003 (5.185) 



Table 2: Results of the fit of the non-perturbative e parameter for the Peterson 
fragmentation function. The value of x 2 /dof is given in parenthesis. The range of the 
fit is indicated in Figs. [T-ll with small crosses. ^We have excluded the first three 
experimental points for the ot 2 s fixed-order fit to the OPAL data. 
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In Figs. [Ty and 11, we illustrate the NLL differential cross section, obtained with 
the exponentiated condition ( |7.1|) . We see a better behavior at x—*l, while the 
values of the e parameter for the best fit are quite small. For comparison, we have 
also plotted the NLL curves at the same e value and the NLL best fit with the value 
of e taken from Tab. |2|. 

The small differences we find in the LL and NLL sectors, with respect to the 
results of Ref. ||, are due to different range, normalization and adjustment of physical 
parameters. 

From the value of the e parameter for ARGUS, OPAL and ALEPH, we can easily 
see that it scales nearly quadratically in the heavy-quark mass, as expected. 



8 Conclusions 

In this work we have considered the heavy-flavour fragmentation functions in e + e~ 
annihilation. We have devised and implemented a method by which all perturbative 
effects that have been calculated so far can be included in the computation of the 
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0.2 0.4 0.6 0.8 1 

X E 

Figure 6: Fragmentation function 1/otot da B /dxE for ALEPH. The value 
of the Peterson e parameter is fixed at 0.0035. 

fragmentation function. These include leading and next-to-leading logarithmic re- 
summation, fixed-order effects up to O (ct^), and mass effects to the same order. 

Our finding can be easily summarized as follows. We generally find little dif- 
ference between our results and the NLL resummed calculation. This indicates that 
mass effects are of limited importance in fragmentation- function physics in e + e~ anni- 
hilation. On the other hand, our calculation confirms the fact that, when NLL effects 
are included, the importance of a non-perturbative initial condition is reduced. For 
example, at ARGUS energies, we see a strong reduction of the e parameter, from 
0.053 to 0.018. This reduction is also observed in the fixed-order calculation, where e 
goes from 0.058 to 0.035 when the 0(a 2 s ) effects are included. 
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Figure 8: Best fit for the improved fragmentation function at OPAL. In 
dashed line, the NLL fragmentation function at the same value of e. 
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Figure 9: Best fit for the improved fragmentation function at ALEPH. In 
dashed line, the NLL fragmentation function at the same value of e. 
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Figure 10: Best fit for the NLL fragmentation function with exponentiated 
starting condition at ARGUS. In dashed line, the NLL fragmentation function 
at the same value of e. In dotted line the NLL best fit. 
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Figure 11: Best fit for the NLL fragmentation function with exponentiated 
starting condition at OPAL. In dashed line, the NLL fragmentation function 
at the same value of e. In dotted line the NLL best fit. 
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